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Abstract 


The  objective  of  the  proposed  research  is  the  investigation  of  pulsed  microwave  and 
combined  microwave/laser  energy  deposition  for  localized  flow  control  in  high  speed  (t.e., 
supersonic  and  hypersonic)  flows.  The  proposed  effort  is  part  of  a  developing  broad  research 
program  in  Electromagnetic  Local  Flow  Control  (ELFC).  The  main  interest  in  ELFC  is 
the  need  for  fast  response  flow  control  systems  for  high  speed  aircraft.  Conventional 
mechanical  or  electro-mechanical  flow  control  devices  have  a  typical  response  time  which 
is  several  orders  of  magnitude  too  large  for  high  speed  aircraft.  ELFC  offers  significant 
potential  for  a  wide  range  of  flow  control  devices  and  techniques  with  fast  response  times. 

There  are  two  overall  objectives  of  the  proposed  research  program.  The  first  objective 
is  to  perform  fundamental  studies  of  microwave  and  combined  microwave/laser  energy  de¬ 
position.  The  studies  include  characterization  of  the  new  microwave  coaxial  resonator  at 
the  University  of  Illinois,  and  developm  ent  of  kinetic-gasdynamic  models  for  microwave 
and  combined  microwave/laser  energy  deposition  in  air.  The  second  objective  is  to  investi¬ 
gate  practical  applications  of  microwave  and  combined  microwave/laser  energy  deposition 
for  flow  control  in  high  speed  flows  relevant  to  the  Air  Force.  These  include  control  of 
1)  supersonic  cavity  flow,  2)  injection  and  mixing  in  a  supersonic  stream,  3)  supersonic 
veliicle  aerodynamics,  and  4)  Edney  IV  interaction. 

The  proposed  research  program  is  an  extension  of  the  current  AFOSR-funded  grant 
entitled  ’'Fundamental  Studies  and  Practical  Applications  of  Electromagnetic  Local  Flow 
Control  in  High  Speed  Flows  Using  Laser  Energy  Deposition”  (AFOSR  Grant  FA9550-04- 
1-0177,  Jan  04  -  Dec  06).  Significant  progress  has  been  achieved  in  the  current  grant  and 
is  summarized  in  the  proposal.  Research  progress  in  the  fundamental  studies  include  de¬ 
velopment  of  a  kinetic  model  for  microwave  energy  deposition  in  air,  a  kinetic-gasdynaniic 
model  for  laser  energy  deposition  in  air,  measurements  of  laser  induced  optical  break¬ 
down.  Research  progress  in  practical  applications  includes  control  of  1)  terminal  (normal) 
shock  waves,  2)  supersonic  flow  over  a  cavity,  3)  Mach  Reflection  to  Regular  Reflection  of 
crossing  shock  waves,  4)  vortex  breakdown  in  shock- vortex  interaction,  and  5)  Edney  IV 
interaction.  In  addition,  two  new  facilities  -  a  microwave  energy  deposition  system  and  a 
supersonic  wind  tunnel  -  have  been  developed  at  the  University  of  Illinois. 

The  proposed  research  program  is  a  collaborative  effort  of  Rutgers  University,  the 
University  of  Illinois  at  Urbana-Chainpaign  and  the  University  of  Minnesota.  In  addition, 
the  research  includes  a  significant  collaboration  with  the  Institute  for  High  Temperatures 
(Moscow,  Russia)  and  St.  Petersburg  State  University  (St.  Petersburg,  Russia)  who  are 
separately  funded  by  EOARD  (London). 
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I.  Introduction 


To  a  large  degree,  the  challenge  of  effective  supersonic  combustion  is  a  problem  of 
mixing.  Aerotherinodynainic  concerns  mandate  a  relatively  short  internal  engine  geometry, 
and  efficient  combustion  requires  fuel  and  air  to  be  well  mixed.  A  fuel  injection  scheme 
must  possess  the  rather  conflicting  qualities  of  low  internal  drag  and  very  rapid  mixing.  The 
problem  is  complicated  by  the  fact  that  flight  experiments  are  extremely  expensive,  and 
generation  of  reasonable  freestream  flow  conditions  in  ground  testing  is  very  challenging. 
In  addition,  accurate  measurement  of  a  liighly  unsteady,  high  speed,  intrinsically  three- 
dimensions  flow  is  difficult.  For  these  reasons,  there  is  great  interest  in  numerical  simulation 
of  supersonic  fuel  injection.  In  order  to  develop  confidence  in  these  numerical  studies,  it  is 
worthwhile  to  consider  even  very  simple  geometries.  One  of  the  most  well-studied  scenarios 
is  the  sonic  transverse  injection  of  fuel  through  holes  flush  with  the  combustor  walls. 

The  gross  flowfield  that  results  from  this  normal  injection  is  well  understood  and 
explained  in  the  literature. 1-7  There  are  many  cases  of  computational  studies;  of  this 
geometry.  Traditionally,  these  have  been  Reynolds-Averaged  Navier-Stokes  (RANS)  simu¬ 
lations  that  result  in  steady-state  solutions.  The  works  of  Chenault  et  al.8  and  Sriram  and 
Matthew9  are  characteristic  of  these  steady-state  nonreacting  RANS  simulations.  These 
studies  are  very  inexpensive  and  capture  some  important  mean  flow  features,  such  as 
the  shock  structures  and  counter-rotating  vortex  pair  (CVP)  seen  in  Figure  1  taken  from 
Reference  10.  However,  these  steady  simulations  are  quite  unable  to  accurately  simulate 
unsteady  turbulent  mixing.  A  more  involved  approach  is  needed  to  capture  these  effects. 

A  hybrid  RANS  and  large-eddy  simulation  (LES)  method  can  be  constructed,  based 
on  the  detached-eddy  simulation  (DES)  formulation  pioneered  by  Spalart  et  al.11  The 
near- wall  region  in  this  approach  Ls  handled  by  the  RANS  portion  of  the  method,  which 
greatly  reduces  the  computational  cost  compared  to  a  LES  resolved  all  the  way  to  the  wall. 
This  wall-modeled  large-eddy  simulation  (WMLES)  approach  allows  routine  computation 
of  flows  at  Reynolds  numbers  beyond  the  reach  of  wall-resolved  LES.  This  methodology 
has  been  previously  applied  to  this  class  of  problem.  Srinivasan  and  Bowersox12  simulated 
several  non-reacting  flows  with  the  unmodified  DES  method.  Boles  et  al.13  investigated 
non-reacting  injection  with  a  hybrid  RANS/LES  method  and  compared  their  results  with 
Mie  scattering  images.  Kawai  and  Lele14  investigated  mixing  mechanisms  of  a  normal 
circular  injector  in  a  LES  simulation.  With  further  modifications  of  the  DES  methodol¬ 
ogy,  Peterson  and  Candler15  investigated  several  reacting  flows  with  the  circular  normal 
injector. 

In  this  report,  we  present  results  from  a  non-reacting  modified  DES  simulation  of  the 
aforementioned  geometry,  run  at  flow  conditions  for  which  Lazar  et  al.16  gathered  particle 
image  velocimetry  (PIV)  data.  This  comparison,  an  extension  of  our  previous  work,  is 
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done  in  order  to  further  validate  the  effectiveness  and  accuracy  of  this  methodology.  The 
simulation  compares  both  mean  and  root-mean-square  (RAIS)  streamwise  and  wall-normal 
velocity  fields  on  the  centerline  at  a  range  of  wall  distances. 

Lazar  et  al.  also  collected  data  on  laser  energy  deposition,  with  the  laser  pulse  located 
at  the  center  of  the  jet  outflow  plane.  High-speed  photography  allowed  the  tracking  and 
investigation  of  the  blast  waves  resulting  from  this  energy  deposition.  We  have  endeavoured 
to  duplicate  this  study  of  pnlse  blast  wave  dynamics.  A  straightforward  and  low-cost 
deposition  technique  has  been  implemented  in  our  simulation  software,  using  a  simple 
local  source  term  in  the  energy  equation.  This  technique  does  not  attempt  to  exactly 
model  the  physics  within  the  volume  of  intense  irradiance.  However  it  is  expected  to 
accurately  model  the  resultant  hot  region  to  the  degree  required  for  high-quality  large- 
scale  simulations.  The  implementation  of  the  energy  deposition  has  not  yet  been  fully 
validated,  and  this  effort  continues  as  of  this  writing. 
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Figure  1.  Illustration  from  Gruber  et  al.10  of  transverse  injection  into  a  supersonic 
crossflow:  (a)  instantaneous  side  view,  (b)  3D  view  of  average  features. 


II.  Numerical  Methods 

The  solver  used  in  this  study  solves  the  compressible  Navier-Stokes  equations  using 
a  hybrid  structured/unstructured,  cell-centered  finite  volume  approach.1""  The  inviscid 
fluxes  can  be  evaluated  using  a  low  dissipation  version  of  the  Steger- Warming  method  with 
either  a  second-order  or  third-order,  upwind-biased  reconstruction.  A  recently  developed 
low-dissipation  flux  evaluation  method  which  is  kinetic  energy  conserving  has  also  been 
incorporated  into  the  flow  solver.  Smooth  regions  of  the  flow  are  handled  by  this  low- 
dissipation  scheme;  in  regions  of  strong  gradients,  the  flux  evaluation  reverts  to  Stegcr- 
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Wanning.  The  gradients  in  the  viscous  terms  are  calculated  using  weighted  least-squares 
fits. 

Time  integration  within  the  solver  is  currently  formally  first-order  accurate.  The 
parallel  Pull-Matrix  Point-Relaxation  method  (FMPR)18  is  used  for  time  integration.  The 
FMPR  method  is  fully  implicit  and  is  not  subject  to  the  unpractical  time  step  limitation 
of  explicit  methods. 

The  unstructured  approach  and  implicit  time  integration  are  enabling  features  for 
simulating  the  type  of  complex  flows  of  interest  in  this  report.  The  unstructured  method¬ 
ology  allows  for  key  flow  features  to  be  well  resolved  without  the  penalty  of  carrying  that 
resolution  to  other  parts  of  the  domain  where  it  is  unnecessary.  That  is,  the  grids  can 
be  made  to  make  every  cell  count.  Implicit  time  integration  allows  for  time  steps  to  be 
taken  which  are  large  enough  that  meaningful  amounts  of  physical  time  can  be  simulated 
in  reasonable  amounts  of  computer  run  time. 

A.  Turbulence  Modeling 

The  simulation  presented  uses  a  wall-modeled  large-eddy  simulation  (WMLES)  ap¬ 
proach  to  turbulence  closure,  with  Reynolds-Averaged  Navier-Stokes  (RANS)  providing 
the  wall  closure.  The  switch  between  wall-modeled  RANS  regions  and  LES  regions  is 
handled  by  a  formulation  developed  by  Travin  et  al.19  The  particular  approach  used  for 
the  simulation  is  an  Improved  Delayed  Detached-Eddy  Simulation  (IDDES)  formulated 
by  the  same  group  and  presented  by  Sliur  et  al.20  The  one-eqnation  turbulence  model  of 
Spalart  and  Allmaras.21  with  density  corrections  of  Catris  and  Aupoix,22  is  used  as  the 
baseline  RANS  model.  The  DES  constant  used  Ls  the  standard  value  of  0.65.  This  IDDES 
formulation  with  the  aforementioned  single-equation  model  shall  henceforth  be  referred  to 
as  the  hybrid-SA  model. 

In  the  interest  of  brevity,  full  details  of  the  hybrid  method  will  not  be  repeated  here;  the 
interested  reader  is  referred  to  the  discussion  by  Travin  et  al.19  However,  several  items  are 
of  note.  DES 97  suffers  from  the  limitation  that  the  LES-RANS  switch  is  solely  a  function 
of  the  grid,  and  so  applies  LES  methodology  in  attached  regions  where  no  turbulence  is 
present.  The  switch  used  in  IDDES  is  a  function  of  the  solution  itself,  so  that  RANS 
behavior  is  maintained  in  regions  where  no  LES  content  is  present.  The  IDDES  method 
also  redefines  the  subgrid  length  scale  near  the  wall,  addressing  the  ‘log-layer  mismatch" 
WMLES  issue  identified  by  Nikitin  et  al.23 

B.  Inflow  Conditions 

The  configuration  of  the  wind  tunnel  nozzle  used  in  the  experiment  is  known,  allowing 
the  aforementioned  upsteain  RANS  calculation  to  generate  the  proper  boundfiry  layer 
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thickness.  The  freestream  stagnation  pressure  in  the  experiment  was  279  kPa  and  the 
stagnation  temperature  300  K;  these  values  were  used  in  the  upstream  calculation.  The 
stagnation  pressure  of  the  injectant  was  378  kPa,  and  the  stagnation  temperature  was  also 
300  K.  The  injector  diameter  was  4.8  mm.  Both  the  injectant  and  the  bulk  flow  were  air. 
The  simulation  was  performed  as  a  two-species  calculation  with  both  species  identical. 

C.  Grid  Generation  and  Boundary  Conditions 

The  grid  used  for  this  simulation  is  composed  of  unstructured  hexahedral  ceLs  (Fig¬ 
ure  2).  It  was  created  using  the  commercial  grid  generation  software  GridPro,  which  is 
amenable  to  sophisticated  topology-based  local  refinement.  This  allows  the  dense  cluster¬ 
ing  of  cells  in  the  region  of  the  jet  plume,  smoothly  and  rapidly  decreasing  in  resolution 
in  three  dimensions,  wit  hout  resulting  in  regions  of  locally  high  anisotropy.  The  computa¬ 
tional  domain  extends  10  injector  diameters  out  from  the  wall,  10  diameters  to  each  side,  5 
diameters  upstream  and  15  diameters  downstream.  The  grid  is  clustered  to  the  wall  such 
that  the  first  cell  was  within  one  wall  unit.  The  final  grid  contains  10,760,218  hexahedral 
elements. 


Figure  2.  Example  of  nested-refinement  grid.  Note  that  grid  is  shown  with  substantial 
global  de-refinement  for  clarity 

The  wall  containing  the  injector,  and  the  injector  walls  and  plenum,  are  adiabatic  no¬ 
slip  walls  with  a  temperature  of  300  K.  The  plenum  contains  a  subsonic  injection  boundary 
condition,  in  which  the  solver  adjusts  pressure  in  order  to  satisfy  a  specified  mass  flow  rate 
and  total  temperature.  The  far-field  walls  are  outflow  conditions,  as  the  region  of  interest 
is  far  from  the  walls  and  showed  no  reflected  shocks  in  the  experimental  results. 
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III.  Results  and  Discussion 

Simulations  were  run  on  256  cores  of  a  quad-core  AMD  Opteron  cluster.  The  time  step 
was  chosen  such  that  the  local  CFL  number  (CFL  number  is  the  non-dimensionalized  time 
step  based  on  the  grid  spacing  and  the  maximum  characteristic  flow  speed)  for  cells  within 
the  jet  plume  was  less  than  unity.  At  the  chosen  time  step,  approximately  174  iterations 
are  required  to  computer  one  flow  time  (based  on  the  jet  exit  diameter  and  the  freestrcam 
velocity).  A  single  iteration  requires  approximately  3.1  seconds  with  the  number  of  cores 
used,  so  that  160  flow  times  can  be  computed  in  a  24  hour  period. 

The  simulation  was  initialized  in  conventional  RANS  mode,  with  conventional  flux 
evaluation  and  large  time  steps,  in  order  to  quickly  establish  the  injector  mass  flow  and 
to  allow  undesired  numerical-startup  transients  to  flow  out  of  the  domain.  This  flow  was 
allowed  to  proceed  for  160  flow  times  (which,  at  these  much  larger  timesteps,  was  completed 
in  a  few  hours)  in  order  to  ensure  that  the  injector  mass  flow  had  reached  steady  state. 
The  simulation  was  theii  restarted  in  DES97  mode,  and  the  low-dissipation  flux  evaluation 
activated.  This  intermediate  step  is  necessary  due  to  the  nature  of  the  hybrid-SA  model, 
which  will  not  switch  regions  to  LES  if  they  do  not  already  possess  LES  content.  Of  course, 
there  is  zero  LES  content  in  pure  RANS  results.  The  intermediate  DES97  iterations  seed 
the  domain  with  LES  content  in  the  appropriate  regions.  This  DES97  simulation  was  run 
for  5000  iterations,  or  about  29  flowrtimes.  The  simulation  was  once  more  restarted,  but 
this  time  in  full  hybrid-SA  mode.  The  simulation  was  then  run  for  50  flowtimes  to  ensure 
that  all  restart-related  transients  have  disappeared,  and  the  simulation  was  run  for  300 
flowrtimes  to  completion.  Statistical  data  were  generated  during  these  last  300  flow  times, 
and  it  is  these  data  that  are  compared  to  experiment. 

A.  Mean  Velocity  Field 

Figure  3  shows  two  time-averaged  streamwise  velocity  fields,  one  from  experiment  and 
one  from  simulation.  The  experimental  images  are  ensemble  averages  of  640  PIV  images; 
the  simulated  velocity  field  is  averaged  over  300  flow  times.  Both  plots  show  some  of  the 
characteristic  flow  details  represented  in  Figure  1,  although  the  upstream  edge  of  the  barrel 
shock  does  not  stand  out.  There  is  only  slight  disagreement  in  velocity  magnitude  in  the 
region  aft  of  the  barrel  shock.  Note  that  the  PIV  data  do  not  extend  precisely  to  the  wall; 
the  method  is  not  amenable  to  measurements  arbitrarily  close  to  a  solid  boundary. 

Figure  4  compares  the  wall-normal  velocity  fields,  averaged  as  previously.  Here,  the 
barrel  shock  and  Mach  disk  are  immediately  apparent.  It  appears  that  the  downstream 
recirculation  region  may  be  furt  her  downstream  in  the  simulation,  but  the  lack  of  near- wall 
PIV  data  does  not  allow  a  concrete  conclusion. 

Overall,  agreement  between  experiment  and  simulation  is  very  good.  The  simulation 
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is  accurately  capturing  all  of  the  flow  physics  that  can  be  resolved  by  the  PIV  data  from 
experiment.  All  of  the  steady-state  centerline  flow  features  identified  in  the  literature,  and 
shown  ill  Figure  1,  are  clearly  visible. 

B.  RMS  Velocity  Field 

In  many  respects,  the  RMS  plots  reinforce  the  observations  derived  from  the  mean 
velocity  fields.  Apparent  in  Figure  5  is  the  difference  in  streauiwise  RMS  velocity  just 
downstream  of  and  below  the  barrel  shock.  The  reason  for  this  disparity  is  not  immediately 
clear,  but  it  should  be  remembered  that  both  experimental  and  simulation  data  jure,  by 
necessity,  spatial  and  temporal  averages.  As  such,  small  differences  in  size  and  location 
of  time-variant  flow  feature  data  does  not  necessarily  indicate  disparity  in  the  actual  flow. 
The  rest  of  the  flowfield  is  in  very  good  agreement.  As  previously  mentioned,  the  large- 
area  variation  very  near  the  wall  is  not  necessarily  noteworthy.  The  areas  of  large'  RMS 
velocities  appear  to  be  somewhat  further  from  the  wall  in  the  simulation,  but  the  near-wall 
limitations  of  PIV  must  be  kept  in  mind. 

Figure  6  does  not  exhibit  any  regions  of  substantial  disagreement.  The  region  down¬ 
stream  of  the  barrel  shock  is  different  in  detail  but  similar  in  general  (especially  keeping 
in  mind  the  space/time  averaging  Issue).  Again,  regions  of  high  RMS  velocities  appear  to 
be  further  from  the  wall  in  the  simulation  compared  to  the  experiment,  particularly  the 
high- valued  red  ‘tail’  trailing  from  below  the  Mach  disk.  The  simulation  also  presents  a 
larger  volume  of  moderate- valued  RMS  velocity  above  the  aforementioned  ‘tail’,  although 
in  both  plots  there  is  a  distinct  lower-valued  region  emanating  from  the  Mach  disk  itself. 
In  general,  the  experimental  and  numerical  data  are  in  good  agreement. 
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Figure  3.  Contours  of  normalized  mean  streaimvise  velocity.  Simulation  upper,  experiment16 
lower. 
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Figure  4.  Contours  of  normalized  mean  wall-normal  velocity.  Simulation  upper,  experiment16 
lower. 
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Figure  5.  Contours  of  normalized  RMS  streamwise  velocity.  Simulation  upper,  experiment 
lower. 
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Figure  6.  Contours  of  normalized  RMS  wall-normal  velocity.  Simulation  upper,  experiment 
lower. 


Simulations  are,  by  nature,  quite  amenable  to  complex  visualization  of  volume  flow 
features.  Figure  7  shows  an  inst  antaneous  plot  of  the  vort  ical  structures  in  the  flow  (demar¬ 
cated  by  the  Q-criterion,  the  second  principle  invariant  of  the  velocity  gradient  tensor).  The 
slices  shown  are  temperature  contours.  Here,  the  flow  features  of  Figure  1  are  extremely 
clear.  Notable  is  the  wide  range  of  length  scales  resolved  by  the  simulation;  capturing 
this  range  of  scales  is  absolutely  critical  to  the  accurate  simulation  of  mixing.  Especially 
prominent  is  the  strong  and  well-defined  horseshoe  vortex  emulating  from  the  upstream 
separation  region.  Also  of  interest  is  the  fact  that  the  inflow  was  not  “seeded”14,25  with 
any  small-scale  fluctuations:  the  highly  complex  upstream  separation  region  is  completely 
a  function  of  the  strong  shock  nearby.  The  low-dissipation  numerics  detailed  above  are  an 
enabling  technology  for  simulation  of  these  structures.  The  physical  generation  and  long 
lifespan  of  these  small  structures  could  not  be  captured  with  other  schemes  without  an 
infeasably  large  grid. 


Figure  7.  Isosurfaces  of  Q-criterion;  temperature  contours. 
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IV.  Conclusions 


Simulations  of  air  injection  into  a  supersonic  crossflow  were  performed  using  a  hybrid 
RANS  /  LES  methodology  based  on  the  detached-eddy  simulation  method,  including  a 
low-dissipation  flux  evaluation  scheme  that  captures  a  large  range  of  length  scales.  This 
methodology  allows  for  relatively  inexpensive  simulation  of  large-scale  unsteady  flow  struc¬ 
tures  in  these  supersonic-combustor-like  geometries.  Results  from  these  simulations  have 
been  compared  to  experiment.  The  mean  velocities  of  the  simulated  flow  were  found  to 
closely  match  the  mean  velocities  from  experiment.  The  RAIS  velocities  exhibited  some 
minor  discrepancies  but  were  very  similar  overall.  Visualization  of  the  entire  volume  of  the 
flow  shows  all  of  the  flow  features  understood  to  exist  in  these  crossflow  injectors.  Laser 
energy  deposition  was  performed  in  experiments.  Numerical  techniques  for  this  energy 
deposition  have  been  developed  and  included  in  the  flow  solver,  but  the  implementation 
has  not  yet  been  fully  validated  and  applied  to  this  case. 


V.  Future  Work 

The  proposed  research  was  to  take  accurate  simulations  of  a  jet  in  crossflow  as  shown 
above,  and  arid  pulsed  energy  deposition  to  model  a  laser  spark.  The  laser  energy  de¬ 
position  local,  energy,  and  frequency  were  to  be  varied  to  sec  if  enhanced  mixing  would 
occur  with  this  means  of  flow  control.  There  have  been  a  number  of  problems  that  have 
prevented  us  from  completing  this  study.  First,  the  experimental  data  were  very  slow  to  be 
produced,  there  were  important  discrepancies  between  simulation  and  experiment  that  had 
to  be  understood  and  solved,  and  the  laser  energy  deposition  model  had  to  be  reconstituted 
and  coded  into  a  different  computational  fluid  dynamics  code.  At  this  point,  all  of  these 
problems  have  been  worked  out  and  we  are  now  beginning  to  produce  meaningful  results 
with  laser  energy  deposition  for  this  problem.  The  results  of  this  work  will  be  presented 
at  the  forthcoming  40th  Amerieal  Institute  of  Aeronautics  and  Astronautics  (AIAA)  Fluid 
Dynamics  conference.26  The  accompanying  paper  will  have  full  details  of  this  study. 
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